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DES EQUATIONS DE DIRAC ET DE SCHRODINGER POUR LA 
TRANSFORMATION DE FOURIER 

JEAN-FRANÇOIS BURNOL 



Résumé. Dyson a associé aux déterminants de Fredholm des noyaux de Di- 
richlet pairs (resp. impairs) une équation de Schrôdinger sur un demi-axe et 
a employé les méthodes du scattering inverse de Gel'fand-Levitan et de Mar- 
chenko, en tandem, pour étudier l'asymptotique de ces déterminants. Nous 
avons proposé suite à notre mise-au-jour de l'opérateur conducteur de cher- 
cher à réaliser la transformation de Fourier elle-même comme un scattering, et 
nous obtenons ici dans ce but deux systèmes de Dirac sur l'axe réel tout entier 
et qui sont associés intrinsèquement, respectivement, aux transformations en 
cosinus et en sinus. 



Dirac and Schrôdinger équations for the Fourier Transform 

Abstract. Dyson has associated with the Fredholm déterminants of the even 
(resp. odd) Dirichlet kernels a Schrôdinger équation on the half-axis and has 
used, in tandem, the Gel'fand-Levitan and Marchenko methods of inverse scat- 
tering theory to study the asymptotics of thèse déterminants. We have pro- 
posed following our unearthing of the conductor operator to seek to realize 
the Fourier transform itself as a scattering, and we obtain here to this end 
two Dirac Systems on the entire real axis which are intrinsically associated, 
respectively, to the cosine and to the sine transforms. 
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1. Abridged English Version 

Dyson jî!| has associated with the Dirichlet kernels 

. sm(s(x-y)) sin(s(x + y)) 

(!) Df(x,y) = r— ± — ■ — r— 

n{x — y) n(x + y) 

on (0, 1) two Schrôdinger équations on (0, +oo) whose potentials are given as: 
(2) W±(s) = -2-^logD±(s)-l 

In this définition the D±(s) are the Fredholm Déterminants àet(ô(x— y) — Df(x, y)) 
associated with the (trace-class) kernels above. After establishing a gênerai link 
between the methods of Gel'fand-Levitan on one hand, of Marchenko on the other 
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hand, and the study of Fredholm déterminants, Dyson obtained in this particular 
case the existence of an asymptotic development as s — > oo of log D±(s), thus 
extending the earlier resuit of des Cloizeaux and Mehta (^J) which is: 

s 2 1 1 

(3) log£> ± (s) = --T glogs + C±+o(l), 

and another heuristic argument enabled Dyson to give the values for the constants 
C±. Except for that last point, still in need of rigorous proof, des Cloizeaux-Mehta's 
and Dyson's results have been later rigorously established by Deift-Its-Zhou ^2] in 
the framework of a Riemann-Hilbert approach (see also Widom ^H] for a proof of 
log£> + (*)£>_( S ) = -4(l+ (l))). 

Our discovery |3 a of the local operator theoretic content of the explicit formulae 
of Number Theory led us to propose to realize the Fourier transform itself as a 
scattering fgj. The multiplicative Fourier analysis of the additive Fourier transform 
leads to two unit modulus functions on the critical line X+( s ) an d X-( s )> the first 
corresponding to the cosine transform, and the second to the sine transform. But 
their phases are of the order 7 log I7I (s = h + l'y) and this escapes the framework 
of scattering "matrices" for potentials which satisfy the usual hypothèses of the 
Gel'fand-Levitan and Marchenko methods. 

We associate in this Note to the cosine transform (resp. sine transform) a dif- 
ferential System which, as in J3j, involves the Fredholm déterminants D±(s). To 
express it as a Dirac System, we shall use, not the variable s = 2na 2 € (0, 00), 
but rather the variable u — log(a) G (—00, +00). Let us thus define D^, a > 
to be the Dirichlet kernel as above for s = 2-ko 2 and let — det(l — D^) and 
d±(u) — da,a = exp(u), u e M. We propose the following Dirac-type System: 



d d 2 
(4a) dû a<yU ^ = ^ydû 2 logd± ^) ' a (") ~ 7/3(«) 



d „, s d 2 



(4b) fcL f3 ^ = + \dû 2 ~ l ° gd± (W) ' ^ + 7 " (U) 



The real number 7 is the spectral parameter. Let /i±(w) = y — -^p logd±(u). 
Each component a(u) and (3{u) has its own Schrôdinger équation. Thèse équations 
are related through a Darboux transformation, and they take the form: 

d 2 

(5a) ~dû? a<yU ^ + ^±( u ) 2 ~ M±( w ))a( u ) = -f 2 a(u) 

(5b) -^ {u) + ( ^ ±(u)2 + ^Wiu) = 7 2 /3(w) 



See also <25at . l|25bll . I|26a|l . Il26b|l . for alternative expressions. The potentials 
arising in l|5âjl . l|5hll are exponentially increasing as u — > +00 and exponentially 
vanishing as u — * —00. We have studied the scattering (rather, reflection) induced 
by <5aH and i|5bjl on waves originating at — 00 and being bounced back to — 00 and we 
have established that the corresponding phase shifts realize (in the manner indicated 
above, and stated more precisely in the French section) the Fourier cosine and sine 
transforms as scattering matrices. A detailed treatment will appear elsewhere. 
Here, we shall explain the origin of the Systems Q and the existence of known 
solutions for each value of the spectral parameter: the structure functions of the 
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Hilbert spaces of entire fonctions which had been associated by de Branges to 
the Fourier transform p]. We gave their first explicit formulae in pÔ] and give 
here their associated differential Systems iQJ: thèse Systems turn out to involve 
the déterminants d±(u), a fact which thus establishes a link between the theory 
of the de Branges Sonine spaces and the quantities arising in the study of the 
random matrices JSj- In [H] we linked through the co-Poisson formula the study 
of the Riemann zeta function to the Sonine spaces. As the reader will perceive, 
référence j2] plays only a distant rôle here. This is despite the fact that we consider 
most beautiful the de Branges discovery [Q of the existence of a natural isometric 
expansion for Fourier self-or-skew-reciprocal functions, and that we are aware of 
the rich contents of his gênerai theory [2] . We are also aware of the rich relations it 
holds with other works, such as the investigations of Krein (see, for example, [Ï3|l. 
relations which are not that easily perceived from [2J. Rather the elucidation of 
the de Branges Fourier Sonine spaces, which was a necessity after 0, was initiated 
in an d is continued here, and turns out to provide an opportunity to partially 
realize some of the ideas which are internai to our approach of the problems of the 
zeta function and which we have expressed elsewhere '.] . [S], jHj- 

2. Un système différentiel 

Nous continuons la discussion entamée dans la partie en langue anglaise du sys- 
tème de Dirac £0 et des équations de Schrôdinger ijÏÏâjl . I|5hjl associées. Nous utili- 
serons les notations suivantes : la transformation de Fourier T agit selon !F(f)(y) = 
J R e 2mxy f(x) dx. On a T = + où T+ est la transformation en cosinus et T- 
est la transformation en sinus. Soit a > 0. Nous utiliserons l'opérateur P a de res- 
triction à l'intervalle [—a, a], ainsi que les opérateurs T^. = J-+P a , = P a JF+P a , 
D+ = (F+) 2 , T~_ = T-P a , F- = P a ^P a , D- = (F-) 2 . Les opérateurs F+ et 
(resp. F a et D a ) peuvent être vus au choix comme agissant sur L 2 (—a,a; ^dx), 
sur son sous-espace pair (resp. impair), ou après un changement d'échelle sur 
L 2 {— 1, 1; icte) ou encore sur L 2 (0, 1; dx). Dans cette dernière représentation le 
noyau de l'opérateur Df est donné par : 

sin(27ra 2 (a; - y)) sin(27ra 2 (x + y)) 

(6) ; \ ± ; t 

•k\x — y) 7r(a; + y) 

Considérons les uniques fonctions ç!>+ et <f)~ solutions (disons, localement inté- 
grables) des équations : 

(7a) 0++^+ = 2cos(27roa:) 

(7b) 0--J^=2cos(27raz) 

Ce sont des fonctions paires, entières. Nous mentionnons également les fonctions 
impaires qui sont associées à la transformation en sinus. Elles vérifient : 

(8a) = 2s,m{2Kax) 

(8b) ^û-^â^â =2sin(27raz) 

Pour être spécifique nous considérerons ici la transformation en cosinus, la dis- 
cussion étant essentiellement identique pour la transformation en sinus. Soit ô x — 
x-§^ + \- Un calcul simple aboutit au système différentiel suivant : 
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(9a) a^-4>+ = 5 x (t> a - (i + a^+{a) + a0 a (a))# 

(9b) a^<t>- = Sxti - (~ - a0+(a) - 

La démonstration consiste àmontrer que aJ^(f)+ —ô x (j)~ et a-ê-(f)~ —àx4>t satisfont 
aussi les équations Q à des facteurs multiplicatifs près qui sont comme indiqués. La 
quantité açi+(a) +a<^>~(a) jouant un rôle important nous la désignerons par /Lt(a) : 

(10) fi{a) = a<j>t (a) + a<j)~ (a) 

On établit sans peine comme conséquence de © : 

et à partir de là la formule suivante : 

(12) M ( a ) 2 = a 2 (0+(a) + C(a)) 2 = a — aj ^ (x)<j>- (x) dx 

3. Le lien avec les déterminants de Fredholm 
Par une formule connue de la théorie générale des déterminants de Fredholm : 

(13) 0+(a) = +^logdet(l+F+) 0"(o) = ~ logdet(l - F+) 

(14) ^+(o) = +^logdet(l + F-) ^-(o) = ~logdet(l-F-) 
Ainsi, on a aussi : 

d , dct(l + F+) 

Nous relions maintenant jjb{a) au déterminant de Fredholm det(l — D+). Par une 
formule bien connue on a : 

(16) ^logd+ = -Tr((l - Dt)- l j- a Dt) 

Or le calcul immédiat du noyau de la dérivée de montre qu'il s'agit d'un opéra- 
teur de rang un, et donc que (1 — D+)~ 1 -£;D+ l'est également. Son vecteur propre 
est donné par la fonction (1 — Z3+) _1 (2 cos(2-7ra 2 a;)) (nous sommes ici sur [0, 1]), et 
on trouve que la valeur propre est égale au produit scalaire de cette fonction avec 
(2a)2cos(27ro 2 x). En faisant usage de la factorisation 1 — = (1 + F+)(l — F+) 
et du caractère auto-adjoint de ces opérateurs on obtient finalement la formule : 

(17) ±i ogd + = - f 4 > +{ x ) ( t>-{x)dx 



où l'on est revenu à l'intervalle [—a, a]. 

La comparaison avec lfÏ2|l . ifïâj) . permet donc d'obtenir le théorème suivant 
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Théorème 1. Soit d+(u) — <i+ = det(l — et H+(u) = fi(a) pour a = exp(u), 

w £ M. On a : 

d 2 

(18) fi+(u) 2 = -—\ogd+(u) 

(19) f i+(u) 2 - f / + (u) = -2^ logdet(l + F+) 

(20) fi+(u) 2 + f i' + (u) = -2^ logdet(l - F+) 



Dans (@J, on a écrit J — J^- logd+(u) pour désigner la fonction u i— ► /J+(w) 



/i(e"). 



4. Le système différentiel des distributions v4 et i? a 



Les fonctions 0+ et 0~ sont aussi des distributions tempérées dont on peut 
prendre la transformée de Fourier. On vérifie de suite que les distributions tempérées 
paires : 

(21a) Aa = ^{T{4> + a ) + 4> + a ) 

(21b) Ba = l -^{Hr a )-<t>a) 

ont la propriété remarquable de s'annuler dans ] — a, a[ tout comme leurs transfor- 
mées de Fourier puisque A a est invariante et B a est anti-invariante. Or un théorème 
général |l()l Thm. 1] nous dit que toute distribution D (paire ou impaire) ayant 
cette propriété a comme transformée de Mellin droite J °° D(t)t~ s dt (a priori défi- 
nie pour Re(s) 0) une fonction entière de s qui a des zéros triviaux en s = —2n, 
n G N (—1 — 2n, dans le cas impair). Nous noterons A a (s) et B a {s) les transfor- 
mées de Mellin de A a et B a , complétées par le facteur Gamma 7r~ s / 2 r(|). Nous 
avons démontré dans ^H] que ces fonctions entières sont en fait les fonctions de 
structure (dans une certaine normalisation) de certains espaces de Hilbert de fonc- 
tions entières que de Branges a associé à la transformation de Fourier QJ. Avec nos 
conventions, l'axe de symétrie de ces espaces est la droite critique, et donc par ce 
fait même, les fonction entières A a et B a vérifient l'hypothèse de Riemann. Nous 
avons aussi montré que la densité de leurs zéros est au premier ordre identique à 
celle pour les zéros de la fonction dzêta de Riemann [Hj. 

Cela étant dit, à partir du système différentiel © on obtient : 

d 

(22a) a da~ Aa = lSxBa ~ ^( a ) A a 

d 

(22b) a ~d~a Ba = ~ l&xAa + ^ a ) B <* 

Au niveau des transformées de Mellin droites on a ô x = s — ^ = ij. On obtient 
donc : 

d 

(23a) a—Aa = ~jB a - fi(a) A a 

oa 

d 

(23b) a— B a = + 7 A + n(a)B a 
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Nous poserons A s (u) — A a (s), B s (u) = B a (s), ^+{u) = fi(a) avec a = exp(u), 
s = i + 17. Ainsi : 

Théorème 2. Les fonctions A s (u) etB s (u) de la variable m € R vérifient le système 
différentiel suivant (s = -z + ij) : 

(24a) ^-A s (u) = - 1 B s (u)- ^ + (u)A s (u) 

du 

(24b) ^-B s (u) = + 1 A s {u) + fi + (u)B s (u) 

du 

Les équations de Schrôdinger associées en découlent : 
(25a) —^A s (u) + ( - 2^ logdet(l + F+ f)A s (u) = 7 2 A s (u) 

(25b) + ( - 2 ^2 ^gdet(l " F a))B s (u) = ~f 2 B s (u) 



Pour la transformation en sinus on obtiendrait 

d 2 , . / d 2 
d^ a{u)+ {- 2 ^ 
d 2 „, , / d 2 



(26a) - — a(u)+ -2 — logdet(l + F") Wu) = 7 2 a(") 



(26b) + ( - 2 ^2 l°gdet(l - *T))/5(«) = 7 2 /3(") 

On notera que si l'on normalise les fonctions de structure A et B des espaces de 
de Branges en imposant A a {\) — 1 alors les équations du deuxième ordre vérifiées 
par A et B prennent la forme d'équations de Sturm-Liouville. 

5. Compléments 

La théorie générale de de Branges [|2] établit l'existence d'équations intégrales 
(qui peuvent parfois prendre, comme il était connu a priori que cela serait le cas 
ici (2|, une forme différentielle) pour certaines chaînes de sous-espaces (dont les élé- 
ments sont aussi des fonctions entières) d'un même espace de Hilbert L 2 (R, dA(p/)). 
Elle les utilise pour obtenir une représentation isométrique de L 2 (R, dA(-j)), ap- 
portant ainsi une vaste généralisation au théorème de Paley- Wiener. Dans le cas 
particulier qui est considéré ici nous avons vérifié que l'expansion isométrique des 
espaces de Sonine (associés à la transformation en cosinus) prend la forme suivante : 
l'application 

f+00 



(27) 



a(u) 
/?(«) 



/(*)= / (A 8 (u)a(u)+B s ( U )0(u)) 



du 



est isométrique de L 2 (M -> M 2 ; 4p) sur L 2 (Re(s) = ^;dA(j)) avec : 
(28) dA( 7 )= , * , a y, s = \+i r 

|7r-*r(|)| 27r 2 

Ce type d'expansion isométrique est aussi intimement lié à des considérations an- 
térieures de Krein, comme il est expliqué dans le livre de Dym et McKean |l.î| . 

Les potentiels pour les équations de Schrôdinger sont exponentiellement crois- 
sants en u — > +00 et s'annulent exponentiellement en u — > —00 : cela suggère, 
au-delà des considérations liées à l'expansion isométrique, d'aborder le problème 
du scattering d'ondes provenant de —00 et y étant renvoyées. Krein a utilisé ses 
techniques pour aborder sur un demi-axe les problèmes de scattering également 
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traités par les approches de Gel'fand-Levitan et Marchenko, dans un cadre où la 
mesure spectrale dA^) satisfait la condition J < oo, ce qui n'est pas le cas ici, 

bien qu'un grand nombre de considérations restent valables. En fait la différence 
essentielle avec la situation habituellement envisagée est qu'ici, heuristiquement, 
plus l'onde est énergétique, plus profondément elle pénétrera dans le potentiel et 
plus longtemps elle mettra pour en être renvoyée : de sorte que l'on s'attend à ce 
que les décalages de phases soient divergents lorsque la fréquence 7 tend vers ±00. 

Ce problème, pour l|25a|l par exemple, est rendu précis par le fait que la condition 
"être de carré intégrable en +00" est une condition aux bords appropriée en +00, 
avec pour unique solution (à une constante multiplicative près) de fréquence réelle 
7 la fonction u 1— > A s {u). Le calcul de la matrice de scattering découle alors de 
la détermination, pour chacune des équations de Schrôdinger l|2f)ajl . j25hjl . Q2fiaf . 
H2fihfl . des solutions satisfaisant la "condition de Jost" en —00. Par exemple, on 
trouve, pour l'équation l|25a,|l . que la solution (complexe) qui se comporte asymp- 
totiquement en —00 comme exp(— iju) est donnée exactement par (a = exp(u), 
s = | + «7, Rc(s) < 1) : 

(29) u ^ e-^ u - e u/2 f tâ{x)x~ s dx 

et à partir de là on établit que la "matrice de scattering" de —00 vers —00 est, pour 
l'équation 125all . exactement donnée par la fonction 7r _ ( 1_s )' 2 r((l— s)/2)/tt~ s / 2 T(s/2) = 
X+(s), s = 9 +*7) Q u i donne la représentation multiplicative (sur la droite critique) 
de la transformation en cosinus. Pour l'équation (|25hl) on trouve — X+{\ + *7) 
comme matrice de scattering. Pour les équations l|26ajl et Il26bll associées à la trans- 
formation en sinus on remplace x+ P ar ia fonction X- donnant la représentation 
multiplicative de la transformation en sinus. Peut-être existe-t-il un moyen d'uti- 
liser cette information pour remonter aux potentiels eux-mêmes et obtenir ainsi 
peut-être une autre approche (spécifique au cas d'un unique intervalle [—a, a]) aux 
résultats de Deift-Its-Zhou et des autres auteurs cités. 

Ces "matrices de scattering" sont aussi les fonctions qui apparaissent dans les 
équations fonctionnelles de la fonction dzêta de Riemann et des séries L de Dirichlet. 
Signalons à ce sujet que la présente Note est la continuation logique de nos travaux 
antérieurs 0, £3 (voir pour un résumé) sur l'opérateur conducteur log |x|+log \y\, 
opérateur qui entre autres choses, révéla que la transformation de Fourier était une 
symétrie des formules explicites. 
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